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In this paper we study theoretically the process of electron capture between one-optical-electron 
atoms (e.g. hydrogenlike or alkali atoms) and ions at low-to- medium impact velocities (v/v e < 
1) working on a modification of an already developed classical Over Barrier Model (OBM) [V. 
Ostrovsky, J. Phys. B: At. Mol. Opt. Phys. 28 3901 (1995)], which allows to give a semianalytical 
formula for the cross sections. The model is discussed and then applied to a number of test cases 
including experimental data as well as data coming from other sophisticated numerical simulations. 
It is found that the accuracy of the model, with the suggested corrections and applied to quite 
different situations, is rather high. 

PACS numbers: 34.70+e, 34.10,+x 



I. INTRODUCTION 

The electron capture process in collisions of slow, highly charged ions with neutral atoms and 
molecules is of great importance not only in basic atomic physics but also in applied fields such as 
fusion plasmas and astrophysics. The process under study can be written as: 

A +q + B~> A 1 -"-^ + B J+ . (1) 

Theoretical models are regularly developed and/or improved to solve (|]J) from first principles for 
a variety of choices of target A and the projectile B, and their predictions are compared with the 
results of ever more refined experiments. 

In principle, one could compute all the quantities of interest by writing the time-dependent 
Schrodinger equation for the system (|l|) and programming a computer to solve it. This task can be 
performed on present-days supercomputers for moderately complicated systems. Notwithstanding 
this, simple approximate models are still valuable: (i) they allow to get analytical estimates which 
are easy to adapt to particular cases; (ii) allow to get physical insight on the features of the problem 
by looking at the analytical formulas; (iii) finally, they can be the only tools available when the 
complexity of the problem overcomes the capabilities of the computers. For this reason new models 
are being still developed 

The present author has presented in a recent paper |J a study attempting to develop a more 
accurate OBM by adding some quantal features. The model so developed was therefore called a 
semi-classical OBM. Its results showed somewhat an improvement with respect to other OBMs, but 
not a dramatic one. 

In this paper we aim to present an OBM for dealing with one of the simplest processes (hi): that 
between an ion and a target provided with a single active electron. Unlike the former one [HI, this 
model is entirely developed within the framework of a classical model, previously studied in pi] (see 
also jjj), but with some important amendments and improvements which, as we shall see, allow a 
quite good accordance with experiments. 

The paper is organized as follows: a first version of the model is presented and discussed in section 
II. In section III we will test our model against a first test case. From the comparison a further 
improvement to the model is proposed (section IV) and tested against the same case, as well as 
other data in section V. It will be shown that predictions with this correction are in much better 
agreement. 
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II. THE MODEL: FIRST PICTURE 



We consider the standard scattering experiment and label T, P, and e respectively the target ion, 
the projectile and the electron. The system T + e is the initial neutral atom. Let r be the electron 
vector relative to T and R the internuclear vector between T and P. In the spirit of classical OBM 
models, all particles are considered as classical objects. 

Let us consider the plane V containing all the three particles and use cylindrical polar coordinates 
(p, z, (j>) to describe the position of the electron within this plane. We can arbitrarily choose to set 
the angle (f> = 0, and assign the z axis to the direction along the internuclear axis. 
The total energy of the electron is (atomic units will be used unless otherwise stated): 

" ^ | - V 2 Z t Z p 

2 2 ^2 + z 2 ^ p 2 + ( R _ z) 2 

Z p and Zt are the effective charge of the projectile and of the target seen by the electron, respectively. 
Notice that we are considering hydrogenlike approximations for both the target and the projectile. 
We assigne an effective charge Zt — 1 to the target and an effective quantum number n to label the 
binding energy of the electron:_E n = Z 2 /2n 2 = l/2n 2 . 

As long as the electron is bound to T, we can also approximate E as 

E(R) = -E n -^ . 

This expression is used throughout all calculations in (I); however, we notice that it is asimptotically 
correct as long as as R — > oo. In the limit of small R, instead, E(R) must converge to a finite limit: 

E(R) -> (Z„ + l) 2 E n 

(united atom limit). For the moment we will assume that R is sufficiently large so that eq . (|^) 
holds, but later we will consider the limit (Q), too. 

On the plane V we can draw a section of the equipotential surface 

U(z,p,R) = -E n -^ . 

This represents the limit of the region classically allowed to the electron. When R — > oo this region 
is divided into two disconnected circles centered around each of the two nuclei. Initial conditions 
determine which of the two regions actually the electron lives in. As R diminishes there can be 
eventually an instant where the two regions become connected. In fig. ^] we give an example for 
this. 

In the spirit of OBMs it is the opening of the equipotential curve between P and T which leads to a 
leakage of electrons from one nucleus to another, and therefore to charge exchange. We make here 
the no-return hypothesis: once crossed the barrier, the electron does not return to the target. It is 
well justified if Z p » 1. As we shall see just below, this hypothesis has important consequences. 

It is easy to solve eq. (|5|) for R by imposing a vanishing width of the opening (p m = 0); further- 
more, by imposing also that there be an unique solution for z in the range < z < R: 

(l + ^T v f - Z v 

K-m — — 



In the region of the opening the potential U has a saddle structure: along the internuclear axis it 
has a maximum at 

z — zq — R- 



<Zp + l 

while this is a minimum along the orthogonal direction. 

Charge exchange occurs provided the electron is able to cross this potential barrier. Let Nn be 
the fraction of trajectories which lead to electron loss at the time t. It is clear from the discussion 
above that it must be function of the solid opening angle angle Q, whose projection on the plane is 
the ±8 m angle. The exact expression for Nq will be given below. Further, be W(t) the probability 
for the electron to be still bound to the target, always at time t. Its rate of change is given by 

dW(t) = -N n dt^-W{t) , 

J- em 



with T em the period of the electron motion along its orbit. 

It is important to discuss the factor dt{2/T em ) since it is an important difference with (I), where 
just half of this value was used. The meaning of this factor is to account for the fraction of electrons 
which, within the time interval [t, t + dt] reach and cross the potential saddle. In (I) it was guessed 
that it should be equal to dt/T em , on the basis of an uniform distribution of the classical phases 
of the electrons. However, let us read again what the rhs of eq. (|^) does mean: it says that the 
probability of loss is given by the total number of available electrons within the loss cone (W(t) x Nn), 
multiplied by the fraction of electrons which reach the potential saddle. However, on the basis of 
the no-return hypothesis, only outgoing electrons can contribute to this term: an electron which 
is within the loss cone and is returning to the target from the projectile is not allowed, it should 
already have been captured and therefore would not be in the set W. It is clear, therefore, that the 
effective period is T em /2, corresponding to the outgoing part of the trajectory. 
A simple integration yields the leakage probability 



Pi = P(+oo) = 1 - W(+oo) = 



+ *71 



= 1-exp (-7^— I N ndt 

\ J- em J — t m 

In order to actually integrate Eq. (^) we need to know the collision trajectory; an unperturbed 
straight line with b impact parameter is assumed: 



R= ^b 2 + (vt) 2 . 

The extrema ±t m in the integral Q are the maximal values of t at which charge exchange can 
occur. If we identify this instant with the birth of the opening, using eq. Q and (^), we find 

VR m - b 2 

— 

V 

At this point it is necessary to give an explicit expression for Nn- To this end, we will consider 
first the case of an electron with zero angular momentum (1 = 0), and then will extend to nonzero 
values. 

In absence of the projectile, the classical electron trajectories, with zero angular momentum, are 
ellipses squeezed onto the target nucleus. We are thus considering an electron moving essentially in 
one dimension. Its hamiltonian can be written as 



P 2 1 



The electron has a turning point at 



-E n 
r 



J_ 

E n 



Obviously the approaching of the projectile modifies these trajectories. However, in order to make 
computations feasible, we make the following hypothesis: electron trajectories are considered as 
essentially unperturbed in the region between the target and the saddle point. The only trajectories 
which are thus allowed to escape are those whose aphelia are directed towards the opening within 
the solid angle whose projection on the V plane is ±# m (see fig. |l|) provided that the turning point of 
the electron is greater than the saddle-point distance: r c > zq. The validity of these approximations 
can be questionable, particularly if we are studying the collision with highly-charged ions, which 
could deeply affect the electron trajectory. We limit to observe that it is necessary in order to make 
analytical calculations. A posteriori, we shall check the amount of error introduced by such an 
approximation . 

The angular integration is now easily done, supposing a uniform distribution for the directions of 
the electrons: 

N n = i(l-cos6> m ) . 

In order to give an expression for 6 m we notice that cos9 m = zo/(p m + zq) 1 ^ 2 , with p m root of 

Z? 2 \ / 7 P 2 

-ft I . r-, I 2 iJptl 



It is easy to recognize that, in the right-hand side, the first term is the potential due to the electron- 
target interaction, and the second is the electron-projectile contribution. Eq. (|l^) cannot be solved 
analytically for p m except for the particular case Z p — 1, for which case: 



2 

Pm = 



/ 2 \ 2 (R s2 



\E(R) 



The form of E(R) function of R cannot be given analytically, even though can be quite easily 
computed numerically ||. In order to deal with expressions amenable to algebraic manipulations, 
we do therefore the approximation: first of all, divide the space in the two regions R <C Ru , R ^ Ru , 
where R u is the internuclear distance at which the energy given by eq. (K|) becomes comparable 
with its united-atom form: 

E n + = i Z V + l) 2 ^n - Ru = ^ ' 



Ru v P ' (Z p + l) 2 -lE n 

We use then for E(R) the united-atom form for R < R u , and the asymptotic form otherwise: 

E(R) = E n + ^, R>R U 
= (Z p + l) 2 E n , R < Ru 
It is worthwhile explicitly rewriting eq. ([li]) for the two cases: 

9%m = R% ( jEjt+W ~ l) ' R>Ru 
' / 1 - R 2 \ R< Ru 



4 \E£ 

and the corresponding expressions for Nn are: 

Na= 1-COBfln, = h 3 _ EnR); R>Ru 
Z o 

= hl-E n R), R<Ru ■ 

Note that Nn = 1/2 for i? = 0. This is a check on the correctness of the model, since, for sym- 
metrical scattering at low velocity and small distances we expect the electrons to be equally shared 
between the two nuclei. 



When Z p > 1 we have to consider two distinct limits: when R — » oo we know that eventually 
p m —> (eq. ^|). It is reasonable therefore to expand ( |l~5| ) in series of powers of p m /R and, retaining 
only terms up to second order: 



2\/ Zp 



-R 2 



— Zp — E„R 



Consistently with the limit R — > oo, we have used the large-7? expression for E(R). 

The limit R — > is quite delicate to deal with: a straightforward solution of eq. (|15|) would give 



Pm 



(Zp + l)E r 



+ 0(R) 



but calculating cos 9 m and eventually Nn from this expression gives wrong results: it is easy to work 
out the result Nn = 1/2, R — > 0. This is wrong because, obviously, the limit Nn -> 1,Z P -> oo 
must hold. The reason of the failure lies in the coupling of eq. ( |l5| ) with the united-atom form for 
E(R): one can notice that the expression thus written is perfectly simmetrical with respect to the 
interchange projectile-target. Because of this symmetry, electrons are forced to be equally shared 
between the two nuclei. This is good when dealing with symmetrical collisions, Z p — Z t = 1, and 
is actually an improvement with respect to (I), where eq. (|2l| ) was used even for small R's and one 
recovered the erroneous value Nn(R = 0) =3/8. But when Z p > 1 the asymmetry must be retained 
in the equations. The only way we have to do this is to extend eq. (|2lj) to small R, obtaining 



It is straightforward to evaluate eq. (|23|) in the limit Z p — > oo, R — > 0, and find the sought result, 2. 
We notice that, from the numerical point of view, it is not a great error using eq. ( [if] ) everywhere: 
the approximation it is based upon breaks down when R is of the order of R u or lesser, which is 
quite a small range with respect to all other lengths involved when Z v > 1, while even for the case 
Zp = 1 it is easy to recover (see equations below) that the relative error thus introduced on Pi is 
APi/Pi = 1/24 for small b (and-obviously-it is exactly null for large b). Therefore, eq. ( pl| ) could 
be used safely in all situations. However, we think that the rigorous altough quite lengthy derivation 
given above was needed since it is not satisfactory working with a model which does not comply 
with the very basic requirements required by the symmetries of the problem at hand. 



We have now to take into account that the maximum escursion for the electron is finite. If we put 
r c — zo and use for zq, r c respectively the expressions given by (Q) and (^), we obtain an equation 
which can be easily solved for R: 



R — R' r 



/ Zp + l)r c 



The R' m thus computed is the maximum internuclear distance at which charge exchange is allowed 
under the present assumptions. Since R'm < Rm (compare the previous result with that of eq. ^| ) 
we have to reduce accordingly the limits in the integration in eq. (^|): it must be performed between 
±t' m , with the definition of t' m the same as t m but for the replacement R m — > R'm- 
The result for the leakage probability is: 



exp 



i -F(Mm) + Gz 

Tern 



where we have defined 



F(u) = 
Gz = 



W Z v + I) 2 



((V^ + i) 2 -^) ~< 



2v 



arcsinh(it] 



(3F(u u ) - 2t u ) (Z p = 1) 
= (Z p >\) , 

Mm = Vt' m /b , 
U u = Vtu/b , 

_ VR U - V 



The period can be easily computed by 

Tern — 2 



dr r- 
— = V2 
V 



l/En 



dr 



= 2-KTi 



En 



(this result could be found also in ||). 

The cross section can be finally obtained after integrating over the impact parameter (this last 
integration must be done numerically): 



2tt 



bPi(b)db 



Again, we have used the fact that the range of interaction is finite: the maximum allowable impact 
parameter b m is set equal to R' m - 



Finally, we consider the case when the angular momentum is different from zero. Now, orbits are 
ellipses whose minor semiaxis has finite length. We can still write the hamiltonian as function of 
just (r,p): 



P_ 

2 



II. 

2r 2 



= —E n 



L is the usual term: L 2 — 1(1 + 1). The turning points are now 



and R' m = + l)r+. 

Now the fraction of trajectories entering the loss cone is much more difficult to estimate. In 
principle, it can still be determined: it is equal to the fraction of ellipses which have intersection 
with the opening. Actual computations can be rather cumbersome. Thus, we use the following 
approximation, which holds for low angular momenta I « n (with n principal quantum number): 
ellipses are approximated as straight lines (as for the I = case), but their turning point is correctly 
estimated using eq. (pfj). Note that also the period is modified: its correct expression is 




III. A TEST CASE 



As a first test case we consider the inelastic scattering Na + + Na(28d, 29s). We investigate this 
sytem since: (i) it has been studied experimentally in j7J; (ii) some numerical simulations using the 
Classical Trajectory Monte Carlo (CTMC) method have also been done on it g|, allowing to have 
detailed informations about the capture probability Pi function of the impact parameter, and not 
simply integrated cross sections; (iii) finally, it has been used as test case in (I), thus allowing to 
assess the relative quality of the fits. 

In fig. (jij) we plot the normalized cross section a = <r/n 4 versus the normalized impact velocity 
v = vn for both collisions nl = 28d and nl = 29s (solid line). The two curves are very close to each 
other, reflecting the fact that the two orbits have very similar properties: the energies of the two 
states differ by a very small amount, and in both cases E n L 2 « 1. The two curves show reversed 
with respect to experiment: er(28d) it is greater than er(29s). The reason is that the parameter r c 
is larger in the former case than in the latter. 

We can distinguish three regions: the first is at reduced velocity around 0.2, where a steep increase 
of cross section appears while going towards lower velocities. Over-barrier models do not appear to 
fully account for this trend: they have a behaviour at low speed which is ruled approximately by 
the 1/v law, consequence of the straight-line impact trajectory approximation: it is well possible 
that this approximation too becomes unadequate in this region. 

The second region covers roughly the range 0.3 -r- 1.0. Here the nl — 29s data are rather well 
simulated while the present model overestimates the data for nl = 28d. The bad agreement for nl = 
28d was already clear to Ostrovsky which attributed it to a deficiency of the model to modelize l- 
changing processes. It seems clear that neither our treatment of the angular momentum is sufficient 
to cure this defect. 

Finally, there is the region at v > I , where again the OBM, as it stands, is not able to correctly 
reproduce the data. The reason for this discrepancy can be traced back to the finite velocity of 
the electron: the classical electron velocity is v e — 1/n, so v can be given the meaning of the ratio 
between the projectile and the electron velocity. When v > 1 the projectile is less effective at 
collecting electrons in its outgoing part of the trajectory (i.e. when it has gone beyond the point of 
closest approach). In simple terms: an electron is slower than the projectile; when it is left behind, 
it cannot any longer reach and cross the potential barrier. 



IV. CORRECTIONS TO THE MODEL 

This picture suggests a straightforward remedy: a term must be inserted in eq. (^|) to account for 
the diminished capture efficiency. This is accomplished formally through rewriting Nn — * w(t, v)Nn, 
with w < 1 . We have put into evidence that w can in principle be function of time and of the impact 
velocity. The simplest correction is made by assuming a perfect efficiency for v < 1, w(t,v < I) = I, 
while, for v > 1, no electrons can be collected after that the distance of minimum approach has 
been reached: w + = w(t > 0, v > 1) = 0. This can appear too strong an assumption, since those 
electrons which are by the same side of the projectile with respect to the nucleus, and which are 
close to their turning point may still be captured. In fig. (Q) we can compare the original data with 
those for w + = (dashed line) . The sharp variation of a at v — 1 is obviously a consequence of the 
crude approximations done choosing w which has a step-like behaviour with v. 

To get further insight, we plot in fig. |^the quantity bPi(b) versus b for the collision Na + +Na(28d). 
The impact velocity is v = 1. The symbols are the CTMC results of ref. |8j. Solid line is the model 



G 



result for w + = 1; dotted line, the result for w + = 0; dashed line, an intermediate situation, with 
w — 1/2. Striking features are, for all curves, the nearly perfect accordance of the value b w 3000 
at which Pi = (it is b m according to our definition). The behaviour at small fe's (Pi m 1/2) is well 
reproduced for w + = 1 while it is slightly underestimated by the two other curves. On the other 
hands, only by setting w + = it is possible to avoid the gross overestimate of Pi near its maximum. 

It is thus evident that the agreement is somewhat improved in the region v ~ 1 by letting 
w + — 0. However, the high-velocity behaviour is still missed by the model, which predicts a power- 
law behaviour a oc while the actual exponent is higher. Within our picture, this suggests 
that also the capture efficiency w~ — w(t < 0) must be a decreasing function of v. An accurate 
modelization of the processes which affect this term is difficult, and we were not able to provide it. 
However, some semi-qualitative arguments can be given. Let us review again the process of capture 
as described in section II and shown in fig. (jl|): if v > 1, an electron at time t can be in the loss 
cone and still not to be lost, since within a time span At w p m /v the position of the loss cone has 
shifted of such an amount that only those electrons which were closer to the saddle point than a 
distance v e At could be caught. The fraction of these electrons is At(2/T em ) ~ p m (2 / vT em ) . This 
correction gives an additional 1/v dependence, thus now a ~ 1/v 2 . 

As an exercise, we try to fit experimental data using to as a free parameter instead that a function 
to be determined by first principles. We choose one of the simplest functional forms: 

1 + \B\ m 

w = ' — ' 

l + \v-j3\ m ' 

with f3, m free parameters to be adjusted. This form gives the two correct limits: w l,v — ► 0, 
and w — > 0, v — > oo. The parameter is not really needed; it has been added to reach a better 
fit. Its meaning is that of a treshold velocity, at which the capture efficiency begins to diminish. 
In fig. (|^) we plot the fit obtained with /3 = 0.2, m = 4 (dotted line): this is not meant to be the 
best fit, just a choice of parameters which gives a very good agreement with data. We see that the 
suggested corrections are still not enough to give the right power-law, if one needs to go to some 
extent beyond the region v = 1. 



V. OTHER COMPARISONS 



A. Iodine - Cesium collisions 



We apply now our model to the process of electron capture 

l q+ + Cs -> l (q - 1)+ + Cs+ 

with q = 6 -T- 30. This scattering process has been studied experimentally in || . It is particularly 
interesting to study in this context since it has revealed untractable by a number of other OBM's, 
including that of (I) (for a discussion and results, see Q). The impact energy is chosen equal to 
1.5 x Z p keV: since it corresponds to v « 1, we can safely assume w = 1. The Cesium atom is in 
its ground state with the optical electron in a s state. 

In fig. ^ we plot the experimental points together with our estimates. In this case the fit is 
excellent. It is important to notice that this agreement is entirely consequence of our choice of 
limiting integration to R given by eq. (|24|): to understand this point, observe that because of the 
very high charge of the projectile, the exponential term in eq. (E5|) is small (F, by direct inspection, 
is increasing with Z p ) and thus Pi w 1. The details of the model which are in F are therefore of no 
relevance. The only surviving parameter, and that which determines a, is R' m . It can be checked by 
directly comparing our fig. Q with fig. 1 of ref. fj, where results from model (I) are shown, which 
differ from ours just in replacing eq. (M) with eq. (0). There, the disagreement is severe. 



B. Ion - Na(n = 3) collisions 

As a final test case we present the results for collisions H-Na(3s,3p). They are part of a set of 
experiments as well as numerical simulations involving also other singly-charged ions: He, Ne, and 
Ar (see Jllj and the references therein and in particular jl2| ; ref. |l^| presents numerical calculations 
for the same system) . In fig. ^| we plot the results of our model together with those of ref. |ll[ . 
Again, we find that only by neglecting w + some accordance is found. The low-energy wing of the 
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curve is strongly underestimated for Na(3s), while the agreement is somewhat better for Na(3p). 
Again, the slope of a for relative velocities higher than 1 could not be reproduced. 
We do not show results for other ions: they can be found in fig. 3 of ref. What is important 

to note is that differencies of a factor two (and even larger for 3s states) appear between light (H + , 
He + ) and heavy (Ne + , Ar + ) ions which our model is unable to predict. We can reasonably conclude 
therefore: (i) that the present model is not satisfactory for v/v e « 1 (it was already pointed out 
in sec. IV) and for v/v e > 1 ; (ii) the structure of the projectile must be incorporated into the 
model otherwise different ions with the same charge should cause the same effect, at odds with 
experiments. As emphasized in [ pd]Jl^ ] the energy defect AE of the process is a crucial parameter: 
captures to states with AE ~ are strongly preferred. Obviously, the value of AE depends on the 
energy levels structure of the recombining ion. 



VI. SUMMARY AND CONCLUSIONS 

We have developed in this paper a classical OBM for single charge exchange between ions and 
atoms. The accuracy of the model has been tested against three cases, with results going from 
moderate-to-good (sec. Ill and IV), excellent (sec. V.A), and poor-to-moderate (sec. V.B). As a 
rule of thumb, the model can be stated to be very well suited for collisions involving highly charged 
ions at low velocities. 

The model is based upon a previous work Q, and adds to it a number of features, which we go 
to recall and discuss: (i) the finite excursion from the nucleus permitted to the electrons; (ii) the 
redefinition of the fraction of lost electrons dt/T em — > dt(2/T em ); (iii) a more accurate treatment of 
the small impact parameter region for symmetrical collisions; (iv) the explicit-altough still somewhat 
approximate-treatment of the capture from / > states; (v) a correction to the capture probability 
due finite impact velocity. Let us discuss briefly each of these points: 

Point (i) and (ii) contribute a major correction: in particular, (i) is essential to recover that excellent 
agreement found in section V.A, while (ii) accounts for the correct bPi behaviour at small fe's (see 
fig. 2). 

Point (iii) is unimportant for actual computations, but corrects an inconsistency of the model. 
Point (iv) has been studied in less detail, in part for the lack of experimental data on which doing 
comparisons. 

Point (v): a good theoretical estimate of w should be of the outmost importance for developing a 
really accurate model of collision at medium-to-high impact velocity. In this paper we have just 
attempted a step towards this direction which, however, has allowed to recover definitely better 
results. 

Finally we recall from sec. V.B that the treatment of the projectile-or better the process of the 
electron-projectile binding-is an aspect which probably awaits for main improvements. We just 
observe that it is a shortcoming of all classical methods, that they cannot easily deal with quantized 
energy levels. 
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FIG. 1. The enveloping curve shows a section of the equipotential surface U = E, i.e. it is the border of the region classically 
accessible to the electron. R is the internuclear distance. The parameter p m is the radius of the opening which joins the 
potential wells, 8 m the opening angle from T; zq is the position of the potential's saddle point. 
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FIG. 2. Cross section for charge exchange for Na + -Na(29s) (upper) and Na + -Na(28d) (lower) collisions. Symbols, experi- 
mental data (adapted from ref. 7); solid line, present model with w + = 1; dashed line, model with w + = 0; dotted line, model 
with w given by eq. (32). Note that the experimental results are not absolutely calibrated, the data shown here are calibrated 
using as reference the CTMC results at v = 1 and nl — 28d. 
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FIG. 3. Probability of electron capture multiplied by impact parameter, Pib, for Na + -Na(28d) collision at v — 1. Squares, 
CTMC data (adapted from ref. 8); solid line, present model with w + = 1; dashed line, w + = 0.5; dotted line, w + — 0. 
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FIG. 4. Cross section for charge exchange in I +9 -Cs collisions. Circles, experimental data with 20% error bar; solid line, 
present model (where we have set w = 1, since we are dealing with v/v e << 1). 
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FIG. 5. Cross section for charge exchange in H + -Na(3s) (upper) and H + -Na(3p) (lower) collisions. Symbols, experimental 
data from ref. (9); lines, present model. 
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